Spin-orbit-driven coherent oscillations in a few-electron quantum dot 
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We propose an experiment to observe coherent oscillations in a single quantum dot with the 
oscillations driven by spin-orbit interaction. This is achieved without spin-polarised leads, and 
relies on changing the strength of the spin-orbit coupling via an applied gate pulse. We derive 
an effective model of this system which is formally equivalent to the Jaynes-Cummings model of 
quantum optics. For parameters relevant to a InGaAs dot, we calculate a Rabi frequency of 2 GHz. 
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Motivated by the desire for a closer understanding of 
quantum coherence and by the drive to develop novel 
quantum computing architecture, a number of break- 
through solid-state experiments have focused on coherent 
oscillations — the back and forth flopping of that most 
fundamental of quantum objects, the two-level system 
HI EJ Q ■ The pioneering work of Nakamura et al. with 
the coherent superposition of charge states of a Cooper- 
pair box Q first demonstrated the possibility of observ- 
ing such oscillations in a wholly solid-state device; thus 
sparking the remarkable progress in qubit development 
in super-conducting systems 0, 13 ■ 

The important double quantum dot experiment of 
Hayashi and co-workers Q showed that coherent oscil- 
lations could also be observed in normal semiconductor 
systems. It is the purpose of this paper to propose an ex- 
periment in which coherent oscillations are observed in 
a single quantum dot (QD), with these oscillations being 
driven by the spin-orbit (SO) interaction. 

The SO interaction in semiconductor heterostructures 
has its origin in the breaking of inversion symmetry, and 
is increasingly coming to be seen as a tool with which to 
manipulate electronic states, see e.g. The grandfather 
of these ideas is the spin-transistor of Datta and Das 
0, in which the SO interaction causes electron spins to 
precess as they move through a two-dimensional electron 
gas (2DEG). In materials where the structural inversion 
asymmetry dominates, e.g. InGaAs, the SO interaction 
can be described by the Rashba Hamiltonian Q 
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In this letter we consider the effects of Hso on electrons 
in a small, few-electron lateral quantum dot. Although 
such dots are yet to be realised in materials with strong 
SO coupling, there is currently a considerable effort to 
develop nanostructures in such materials P . Our interest 
here is not in open or chaotic QDs 0, llfj j. but rather in 
small dots in the Coulomb blockade regime. 

Such dots have been studied by a number of authors 
jlll IT3. IfH|. but our analysis differs in a crucial respect: 
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by making an analogy with quantum optics, we are able 
to derive an approximate Hamiltonian that captures the 
essential physics of the dot. This model is formally iden- 
tical to the Jaynes-Cummings ( JC) model 0| > first de- 
rived in the context of the atom-light interaction. Here, 
the roles of the atomic pseudo-spin and light field arc 
played by the spin and orbital angular momentum of 
the electron respectively. The system then naturally de- 
composes into a set of two- level systems (TLS), any of 
which may be considered as the qubit degree of free- 
dom within which coherent oscillations can occur. These 
oscillations are genuine Rabi oscillations ^5}; with or- 
bital and spin degrees of freedom exchanging excitation. 
This "spin-orbit pendulum" behaviour has been noted in 
three-dimensional models in nuclear physics [lfij . 

Having elucidated the origin and properties of the TLS, 
we then describe an experimental scheme through which 
the coherent oscillations can be investigated. The key 
problem here is that of injecting into, and reading out 
from, states which arc not cigenstates of the SO cou- 
pled system. In the Hayashi experiment Q, this was 
achieved through the spatial separation of the two dots, 
which makes the leads couple to the localised left and 
right states, rather than to the bonding and anti-bonding 
eigenstates. In our single dot system, the direct anal- 
ogy of this would be the injection of spin polarised elec- 
trons. Given the difficulty of interfacing ferromagnetic 
leads with semiconductors Q , we avoid their use by ex- 
ploiting the fact that the strength of the SO interaction 
can be controlled by external gates 0, 0, 0] . 

Our starting point is the Fock-Darwin theory of a sin- 
gle electron in a 2DEG with parabolic confinement of 
energy hu> H3], 
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where m is the effective mass of the electron. Applying a 
perpendicular magnetic field in the symmetric gauge, in 
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second quantised notation we have 



Hn 



with 
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l>q + u> 2 /4. Introduction of 
a± = 2~ 1 / 2 (a x =F ia y ) decouples the system into eigen- 
modes of frequency cu± = ui ± uj c /2. 

We now include the Rashba interaction of Eq. , for 
which the coupling strength a is related to the spin pre- 
cession length /so = h 2 /2ma. With magnetic length 
Ib = ^Jh/mujc, we have 



so 



7 + (a + cr + + a^(T_) - 7_(a_c_ 
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with coefficients j± = 1 ± ^ (J/'sJ and Z = ^h/muj. 

Adding the Zccman term, in which we take g to be neg- 
ative as in InGaAs, performing a unitary rotation of the 
spin such that a z — > — a z and ct± — > — c^, and rescaling 
energies by fkoo we arrive at the Hamiltonian 

H = uj + a < + a + + ui-aLa,- + —E z a z 
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7_(a_cr + + a}_<7_) - 7+(a + er_ + a} + a + ) ,(5) 



where = \A/mwo is the confinement length of the dot 
and E z = \g\m/(2m c )(lB/lo) 2 is the Zeeman energy with 
m e the bare mass of the electron. 

This single-particle picture is motivated by the good 
agreement between Fock-Darwin theory and experiment 
in the non-SO case [2(| , and by studies which have shown 
that many-body effects in QDs play only a small role at 
the magnetic fields we consider here |ll|, |l2j, |2l| . 

We now derive an approximate form of this Hamilto- 
nian by borrowing the observation from quantum optics 
that the terms preceded by 7+ in Eq. (|SJ) are counter- 
rotating, and thus negligible under the rotating-wave ap- 
proximation |l5l | when the SO coupling is small compared 
to the confinement. This decouples the uj + mode from 
the rest of the system, giving H = u) + n + + Hjc where 

Hjc(a) = tu-(J_a- + -E z a z + A(a_(T + + aL<7-), (6) 

with A = Iq^-/211so- This is the well-known Jaynes- 
Cummings model (JCM) of quantum optics. It is com- 
pletely integrable, and has ground state \0,[) with en- 
ergy Eq = —E./2 independent of coupling. The rest of 
the JCM Hilbcrt space decomposes into two-dimensional 
subspaces {\n, t), \n + 1, J.}; n = 0, 1, . . .}. Diagonali- 
sation in each subspace gives the energies E^ 1 '^ = 
(n + 1/2) ui- ± A n /2 with detuning <5 = u>- — E z and 
A„ = yj8 2 + 4A 2 (n + 1). The eigenstates are 
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FIG. 1: Spectral features of Rashba-coupled quantum dot as 
function of magnetic field. The parameters used are typical 
of InGaAs: g = —4, m/m e — 0.05 with dot size Iq = 150 nm. 
Resonance occurs at Bo = 90 mT. (a) Low-lying excitation 
spectrum for spin-orbit coupling a = 0.8 x 10 -12 eVm. (b) 
Lowest lying anticrossing. Thick line is JC model showing 
anticrossing width Ao at 5 = 0, and thin line is exact nu- 
merical result, (c) Plot of width A„ against central energy 
of anticrossing with the dot on resonance for different a in 
the range 0.3 — 2.0 x 10~ 12 eVm. The exact numerical re- 
sults (circles) show excellent agreement with the square-root 
behaviour predicted by the JC model in this a range. 



with tan6»i' l,±) = (6 ± A„)/2A\ATTT. 

Figure ^ shows a portion of the excitation spectrum 
obtained by exact numerical diagonalisation for a typical 
dot in InGaAs. The approximate Hjc describes the en- 
ergy levels of the system to within 10% of the typical an- 
ticrossing width and 1% of luq. This small discrepancy is 
visible in Fig. ^5. In the following, we are only concerned 
with the lowest-lying energy states in the dots. Without 
SO interaction, these states are described by n + = 
indicating that the states converge to the lowest Landau 
level in the high-field limit, and by n,_ corresponding to 
the quantum number of angular momentum. The SO 
interaction thus couples two states of adjacent angular 
momentum and opposite spin. The detuning 8 uniquely 
identifies ui c for fixed material parameters and dot size. 

Under the assumptions of the constant interaction 
model |20l | , the most important prediction of this model 
for linear transport is that, with the dot on resonance, the 
addition-energy spectrum for the first few electrons (up 
to 18 here) is described by a sequence of well-separated 
anticrossings, the width of which increases as a\/n + 1. 
This behaviour is shown in Fig. and its observation 
would be confirmation of our JC model, and would per- 
mit a determination of a in quantum dots. 

We now describe the procedure for observing spin-orbit 
driven Rabi oscillations. Our proposal is somewhat sim- 
ilar to that of Nakamura with a voltage pulse driving 
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FIG. 2: Configuration of the dot in the various stages of 
the cycle, (a) The positions of the dot levels ip^-i , chemical 
potentials [il,r, and the tunnelling rates Fl > Fr. (b) The 
coupling is initially asi. On average, for times ij > F~^ the 
dot will be initialised in state ipai- ( c ) The applied voltage 
pulse lowers the dot levels and non-adiabatically changes the 
coupling to «2 7^ on, thus inducing Rabi oscillations, (d) 
Pulse is switched off after time t p and the levels return to their 
initial places. Tunnelling to right occurs when the electron 
has oscillated into upper state. Relaxation rates Tig are also 
shown. 



the system, but with the crucial difference that the oscil- 
lations here are induced, not by a change in the detuning, 
but by a change in the SO coupling strength. We operate 
in the non-linear transport regime and address a single 
two-level system by being near resonance and by tuning 
the chemical potentials of the leads close to the n-th an- 
ticrossing. The SO coupling is set to ai and the states 
taking part in the oscillation are eigenstates of -ffjc(ai), 
namely ip^ , which are situated symmetrically around the 
chemical potential of the right lead [Ir, see Fig. 01. The 
temperature is taken smaller than the detuning k^,T <C 5 
to avoid the effects of thermal broadening. Assuming 
Coulomb blockade and considering first-order sequential 
tunnelling only, electrons can either tunnel from the left 
lead into the dot via state tp* and subsequently leave to 
the right or, alternatively, tunnel to state tp~ blockad- 
ing the dot, see Fig. Eh- Assuming tunnelling through the 
left /right barrier at a constant rate T^/r, we set > T r 
to assure that the dot is preferentially filled from the left; 
thus maximising the current. On average then, the dot 
will be initialised in state ip~ for times t{ > . 

Having trapped an electron in this state, we apply 
a voltage pulse to the gate. This has two effects. 
Firstly, this change in voltage alters the SO coupling to 
a new value oli. Since this change is performed non- 
adiabatically, the electron remains in the initial eigen- 
state i/>~ until Rabi oscillations begin between this state 
and ip^ under the influence of the new Hamiltonian 



Hjc{ot2)- Secondly, the TLS is drawn below both chem- 
ical potentials, assuring that oscillations can occur with- 
out tunnelling out of the dot, see Fig.Et. 

After a time t p , the gate voltage is returned to its initial 
value, and the TLS resumes both to its original position 
and coupling a±, as in Fig. [2Ji- Tunnelling out of the 
dot can now occur, provided that the electron is found in 
the upper state, which happens with a probability given 
by the overlap of the oscillating wave function at time t p 
with the upper level, 

P(t P ) = l«l*(ip))| 2 = KV'ile-^^^-)! 2 . (8) 

This process is operated as a cycle and the current 
is measured. From probability arguments we see that 
/ as eTnP{tp), where we have used the simplification 
that r^, 1 > t p , r^ 1 . Thus, by sweeping t p we are able to 
image the time evolution of Rabi oscillations, just as in 
the previous experiments of Nakamura and Hayashi. 

The singular case of a non-adiabatic change in a from 
zero to a finite value produces oscillations with the maxi- 
mum possible amplitude, P max = 1- However, in realistic 
systems only changes between finite values of a are fea- 
sible. This leads to a reduction in the amplitude, and 
achieving a significant oscillation signal requires a suit- 
ably large change in a. In experiments with 2DEGs, 
changes in a of a factor of 2 are reported, and in a recent 
Letter by Koga et al, a was shown to vary in the range 
« (0.3 - 1.5) x 10~ 12 eVm (a factor of 5) in one InGaAs 
sample 0| . Grundler 0] has shown that the large back- 
gate voltages usually used to change a can be drastically 
reduced by placing the gates closer to the 2DEG. Thus, 
it is conceivable that changes in a of a factor between 
2 and 5 could be produced with voltages small enough 
to be pulsed with rise times substantially shorter than a 
typical coherent oscillation period. 

In Fig. we plot time-traces of the transition prob- 
ability P(t p ) calculated for the first anticrossing as a 
function of magnetic field. We have used the values 
ai = 1.5x 10 -12 and ct2 = 0.3x 10 -12 eVm from the Koga 
experiment fl9| . The amplitude of the oscillations P max 
for different ratios of ct^/cti is presented in Fig.[3|3, which 
shows a node at B = B$ (6 = 0). This is because, for 
5 = 0, the eigenstates of JCM are 2" 1 / 2 (|n, f)±|n+l, |)) 
for all a 0. Therefore, a finite detuning is required 
to obtain the maximum amplitude, which concurs with 
5 > fceT, Tr to overcome broadening effects. Both the 
amplitude P max and frequency f2 show non-trivial de- 
pendencies on <x\ and ai as well as on the magnetic field. 
This latter behaviour stems from the parametric depen- 
dence on B of all three parameters in -ffje- 

For our model parameters with o^/ai = 1/5 and 
with the detuning set such that the amplitude is max- 
imised, we have P ma x ~ 0.45 with a Rabi frequency of 
$1 = 2 GHz, which corresponds to a period of about 3 ns. 
This is within accessible range of state-of-the-art exper- 
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FIG. 3: Characteristics of the Rabi oscillation, (a) Prob- 
ability P(t p ) of finding electron in upper level after time 
t p following the non-adiabatic change ai = 1.5 — > ct2 — 
0.3 x 10~ 12 eVm as function of magnetic field, (b) Ampli- 
tude of oscillation as function of B/Bo for changing from 
ai = 1.5, 0.8, 0.6 to a 2 = 0.3 x 10~ 12 eVm (top to bottom), 
(c) Phonon-induced relaxation rate for InAs parameters a — 
1.5 x 10" 12 eVm, P = 3.0 x 10~ 21 J 2 /m 2 , p = 5.7 x 10 3 kg/m 3 , 
c = 3.8 x 10 3 m/s. Close to Bo the rate is suppressed to 

r ep < io~ 7 o;o- 



imcntal technique. Note that the period can be extended 
by using weaker confinement and SO coupling. 

For both the observation of coherent oscillations, and 
the operation as a qubit, it is essential that the lifetime 
of state ifj^ i s long. This is the case for a pure elec- 
tronic spin in a QD 0, , and we now show that the 
hybridisation of the spin with the orbitals, and the ensu- 
ing interaction phonons, does not affect this. We assume 
a piezo-electric coupling to acoustic phonons via the po- 
tential V ep = A q e lq r (6 q + b_„), with phonon operators 
| Aq | 2 = hP/2pcqV, with coupl ing P, mass density 



6 q and 

p, speed of sound c, and volume V |23| 
Golden Rule calculation yields the rate 
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For n = 0, 



a 5 m, (9) 



with lu s = c/l , £ = 2- 1 / 2 (l/l )(A/huj s ), and /(£) < 
8/15. Close to B , £ <C 1, and thus the rate is extremely 
small r ep ~ 10 4 s _1 (Fig. Therefore, the robustness 
of spin qubits is not significantly weakened by the SO 
hybridisation. 

In general, residual relaxation affects our measurement 
scheme in two ways. During the oscillation (Fig. Eh), the 
system may relax to the eigenstate ip~ . This damps 
the oscillation by a factor exp(— Tit p ) to the constant 
value / = eTnP mlix /2. Relaxation during the read-out 
phase (Fig. simply reduces the overall amplitude of 
the signal by a factor exp(— r2/r^). Clearly then, to 
observe oscillations, we require Ti < and T 2 < r^. 

In summary, we have outlined a proposal for the ob- 
servation of spin-orbit driven coherent oscillations in a 



single quantum dot. We have derived an approximate 
model, inspired by quantum optics, that shows the oscil- 
lating degree of freedom to represent a novel, composite 
spin-angular momentum qubit. 
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